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1.  INTRODUCTION 


Despite  perennial  research  in  pressure-atomized  liquid  jets,  the  ability  to  characterize  the  core  structure 
of  such  jets  has  been  extremely  limited.  This  limitation  arises  from  the  inability  to  observe  such  structures 
using  standard  visualization  techniques.  For  "atomizing"  jets  (i.e.,  jets  that  begin  to  diverge  immediately 
at  the  nozzle  exit),  the  liquid  core  is  shrouded  by  an  optically  dense  cloud  of  droplets.  Two  experimental 
techniques  have  been  developed  to  overcome  this  obstacle:  an  electric  grid  method  (Chehroudi  et  al. 
1985)  and  flash  x-ray  photography  (Baev  et  al.  1986;  Waiken  and  Krehl  1990).  The  former  identifies  the 
length  of  the  core  penetration  when  the  core  creates  an  electrical  path  between  the  nozzle  and  a 
downstream  electrical  probe.  Despite  the  success  reported  in  obtaining  data  with  this  technique,  it  does 
not  provide  the  overall  core  structure  and  the  measurement  is  intrusive.  Flash  x-ray  photography  takes 
advantage  of  the  fact  that  x-rays,  defined  as  radiation  with  wavelengths  in  the  range  0. 1  to  10  nm,  are  only 
scattered  at  the  molecular  level  following  absorption,  not  at  interfaces.  Thus,  radiography  provides  a 
shadowgraph  corresponding  to  molecular  density.  Despite  the  promise  of  this  technique  and  its  advantages 
over  the  electric  grid  method,  only  a  few  studies  have  been  reported  and  these  studies  have  dealt  primarily 
with  circular  jets  under  nonevaporating  conditions. 

This  report  describes  the  implementation  of  the  flash  x-ray  technique  for  studying  the  core  structure 
of  atomizing  jets  at  high  temperatures  and  pressures.  This  work  stems  from  our  efforts  to  resolve  issues 
related  to  the  development  of  liquid  propellant  gun  technology  (Morrison,  Knapton,  and  Bulman  1988). 
The  current  technology  is  based  on  regeneratively  injecting  a  liquid  monopropellant  into  the  combustion 
chamber.  The  combustion  intensities  in  such  devices  are  extremely  high  and  large-scale  pressure 
fluctuations  have  been  observed  in  nearly  all  such  gun  designs  tested  to  date.  The  pressure  fluctuations 
are  presumably  associated  with  the  unstable  combustion  of  unreacted  liquid  in  the  chamber.  Though  there 
is  no  direct  experimental  evidence,  it  is  considered  that  little  liquid  exists  beyond  the  jet’s  core.  If  this 
is  the  case,  then  the  core  represents  the  bulk  of  the  unreacted  liquid  and  its  distribution  in  the  chamber. 
Thus,  there  may  be  a  correlation  between  the  core  structure  (dynamics)  and  the  observed  pressure 
fluctuations.  Beyond  this  application,  the  current  study  has  ramifications  for  diesel  and  liquid  rocket 
engines,  both  of  which  involve  injection  into  high  temperature  and  pressure  environments.  In  recent  years, 
interest  in  improving  the  performance  of  such  systems  has  led  to  designs  in  which  the  liquid  traverses  its 
critical  point  prior  to  combustion.  There  is  a  tangible  need  for  better  data  on  such  systems.  For  example, 
theory  and  experimental  data  obtained  with  the  electric  grid  method  suggests  the  core  penetration  distance 
(L(^)  can  be  predicted  via  the  equation. 
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Lc  =  CoCIoIpl/Pg)  '  > 


(1) 


where  is  a  constant,  is  the  nozzle  diameter,  and  Pl  and  Pq  are  the  liquid  and  gas  densities, 
respectively.  However,  numerical  simulations  of  diesel  engine  spray  combustion  underestimate  the 
experimentally  observed  spray  penetration  distance  (Gonzalez  and  Reitz  1991). 

The  implementation  of  flash  radiography  for  characterizing  the  core  structure  of  atomizing  jets 
evaporating  at  high  temperatures  and  pressures  is  challenging  because  the  injection  must  be  radiographed 
through  chamber  windows  that  will  not  fail  under  such  conditions.  The  chemical  composition  and 
thickness  of  windows  having  the  required  physical  strength  can  significantly  attenuate  the  x-ray  beam 
intensity.  Also,  the  injection  orifice  must  be  relatively  small  so  that  the  entire  core  can  be  observed  in 
a  reasonably  sized  pressure  vessel.  Both  of  these  requirements  limit  the  contrast  achievable  in  a  single 
flash  x-ray  image.  We  have  overcome  these  difficulties  through  judicious  exp)erimental  design.  Reported 
here  are  the  results  obtained  from  experiments  where  methyl  iodide  (CH3I)  was  injected  at  approximately 
75  m/s  through  a  circular  orifice  into  either  (hot)  nitrogen  or  the  post-combustion  products  of 
H2/Air/Argon  mixtures.  The  results  from  similar  tests  involving  CH3I  injections  at  40  m/s  through  annular 
orifices  are  also  presented.  Methyl  iodide,  an  unusual  choice  for  such  a  study,  has  a  combination  of 
desirable  properties:  namely,  a  high  x-ray  absorption  cross  section,  well-characterized  physical  parameters, 
and  a  critical  point  which  can  be  exceeded  in  our  test  chamber.  In  addition  to  flash  x-ray  data,  jet  images 
based  on  visible  light  scattering  were  obtained  for  injections  into  nitrogen.  The  results  are  analyzed  in 
terms  of  variation  in  gas  density,  temperature,  and  pressure. 

2.  DESIGN  CONSIDERATIONS  FOR  FLASH  X-RAY  EXPERIMENTS 

The  two  most  important  parameters  controlling  the  ability  to  resolve  details  in  a  flash  x-ray  image  are 
the  contrast  and  penumbra.  The  "intrinsic"  contrast  (C)  is  defined  in  terms  of  differences  in  transmitted 
beam  intensity  for  different  thicknesses  of  the  same  material.  The  intrinsic  contrast  is  rendered 
experimentally  by  exposing  a  detector.  (In  this  work,  the  detector  is  a  film-intensifying  screen 
combination.)  This  produces  a  monochrome  image  (transparency)  whose  optical  densities  correspond  to 
transmitted  intensities.  Penumbra  is  an  intrinsic  "unsharpness"  due  to  the  fact  that  the  x-ray  beam  diverges 
from  the  source  (i.e.,  the  beam’s  rays  are  not  parallel  as  they  pass  through  the  object).  Each  of  these 
considerations  is  discussed  in  the  following  sections.  This  discussion  is  based  on  Krehl  and  Warken 
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(1990)  and  a  Hewlett-Packard  Technical  Bulletin  (1973)  and  they  are  recommended  for  readers  interested 
in  a  more  complete  discussion  of  these  aspects  of  radiography. 

For  a  monochromatic  beam  with  wavelength  X,  passing  through  a  sample  of  thickness  x,  the  ratio  of 
transmitted  GJ^l)  to  incident  [IjfX)]  intensity  is  given  by 

T(X)  =  It(X)/Ii(X)  =  exp[-p(X)x]  ,  (2) 

where  p(X)  is  the  linear  absorption  cross  section  of  the  sample,  and  Ij(X)  is  equal  to 


W  =[nTn(^)lls(^)  • 

n 


(3) 


where  l5(X)  is  the  output  of  the  x-ray  source  and  T„(X)  is  the  transmissivity  through  non-object  elements 
in  the  path.  Neglecting,  for  the  moment,  film  response  vs.  exposure  considerations,  it  can  be  shown  that 
the  contrast  [C(X)]  between  two  objects  (n  and  n  -h  1)  of  the  same  material  is  given  by 

C(X)  =  [l„(X)-I„^i(X)/I„(X)]  =  l-exp[-p(X)Ax]  ,  (4) 

where  Ax  is  the  difference  in  thickness  of  objects  n  and  n  +  1.  The  transmission  properties  of  uniformly 
thick,  nonobject  elements  (such  as  windows)  are  not  important  at  this  level  of  analysis.  For  a 
polychromatic  source,  the  contrast  will  be  a  weighted  average  of  the  intensity  distribution, 

C  =  l-/Ii(X)exp[-p(X)Ax]  dX  /  /Ii(X)  dX  .  (5) 

Thus,  designing  an  experiment  to  optimize  the  intrinsic  contrast  requires  a  knowledge  of  the  spectral 
output  of  the  x-ray  source,  the  transmissivity  characteristics  of  the  object  and  nonobject  elements,  and  the 
wavelength  dependence  of  the  sample’s  absorption  cross  section. 
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The  transmissivity  of  nonobject  elements  in  the  path  is  also  important  when  film  response  is 
considered.  For  a  given  film,  optimal  rendering  of  the  intrinsic  contrast  depends  primarily  on  maximizing 
the  change  in  the  optical  density  [OD  =  log  (1^  / 1^)]  per  change  in  x-ray  intensity  incident  on  the  film  (l)\ 
i.e.,  dOD(I)  /  dl.  (Note,  that  in  the  definition  of  OD,  the  parameters  and  are  incident  and  transmitted 
intensities  associated  with  scanning  the  transparency.  They  should  not  be  confused  with  incident  and 
transmitted  x-ray  intensities  associated  with  exposing  the  film.)  In  general,  there  win  be  an  optimal 
incident  x-ray  intensity  (exposure)  range  which  maximizes  dOD(I)  /  dl. 

The  penumbra  (y)  is  a  function  of  the  object-to-film  distance  (fof)*  x-ray  tube-to-film  distance  (r^f), 

and  the  focal  diameter  of  the  x-ray  source  (D^, 

y  =  Dx  r^f  /  (rsf-r^f)  .  (6) 

The  focal  diameter  is  fixed  by  the  source  hardware,  but  penumbra  can  be  minimized  by  decreasing  r^f  and 
increasing  r^^ 

Equations  3,  4,  and  5,  coupled  with  film  response  and  penumbra  considerations,  guided  our  effort  to 
design  a  system  where  the  details  of  the  liquid  core  stracture  could  be  resolved  down  to  0.1  mm.  The 
flash  x-ray  unit  available  to  us  was  a  standard  150-KV  Hewlett-Packard  system.  The  flash  duration  of 
this  system  is  sufficienfly  short  (50  to  70  ns)  to  freeze  the  jet  motion.  Several  different  x-ray 
targetAvindow  (tube  material)  configurations  are  available  for  this  system  and  these  configurations 
determine  l5(X).  Since  we  were  interested  in  resolving  the  details  of  "thin"  objects  and  x-ray  absorption 
cross  sections  generally  increase  with  decreasing  photon  energy — a  spectrum  of  "soft"  (20  keV)  to 
"medium"  (20  to  50  keV)  x-ray  ouq)ut  was  desired.  This  dictated  the  use  of  a  beryllium  (vs.  aluminum) 
windowed  tube  and  a  tungsten  target.  The  tungsten  target  has  an  emission  spectrum  heavily  weighted  with 
"very  soft"  (<  10  keV)  and  medium  x-rays.  In  choosing  test  chamber  windows  and  a  liquid  to  inject, 
absorption  cross  section  was  the  primary  consideration.  Since  absorption  cross  sections  in  the  x-ray 
energy  range  tend  to  increase  with  increasing  atomic  number,  we  searched  for  (1)  liquids  with  a  high 
density  of  high  (>  50)  atomic  number  constituents  and  (2)  polymers  with  low  atomic  number  constituents 
for  windows.  The  search  for  an  appropriate  liquid  centered  on  iodine  (atomic  number  =  53)  containing 
molecules.  CH3I  was  selected  because  it  had  the  desirable  properties  of  being  a  liquid  at  room  conditions, 
a  thermodynamic  critical  point  which  could  be  exceeded  in  our  test  chamber,  and  well-characterized 
physical  properties  (see  Table  1).  Polycarbonate  (Lexan'"''^,  [CjgHj403)„]  and  polypropylene  [(C3Hg)j,] 
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Table  1.  Test  Liquid  (CH3I)  Properties® 


Pl, 

(g/cm^) 

Pl 

(g/cm-s) 

% 

(g/s^) 

P  . 

^cnt 

(MPa) 

'^crit 

(K) 

Oh** 

2.279 

0.005 

25 

7.366 

528 

0.002 

®  From  Handbook  of  Giemistrv  and  Physics.  (Near  room  temperature). 
^  Ohnsorge  number  =  (given  for  dQ  =  0.1  cm). 


were  selected  as  window  materials  based  on  their  ready  availability  and  previous  experience  with  such 
windows  at  the  pressures  of  these  experiments.  These  choices  determined  the  basic  experimental 
conflguration,  and  all  other  considerations  were  based  on  achieving  the  best  results  with  this  hardware  and 
liquid. 

To  render  the  x-ray  image,  a  number  of  film-intensifying  screen  combinations  were  tested.  The 
combination  which  yielded  the  highest  OD/exposure  for  the  basic  experimental  configuration  was  NDT  57 
(Dupont)  film  with  a  Trimax  8  (3M)  intensifying  screen.  For  this  combination,  background  optical 
densities  of  approximately  3.5  were  achievable.  (Note  that  background  areas  correspond  to  the  highest 
optical  densities  in  the  image.)  Such  exposures  would  have  maximized  the  dynamic  range  of  gray  levels 
in  the  image.  Unfortunately,  our  optical  scanner  was  unable  to  reliably  digitize  transparencies  with 
OD  >  2.5.  Therefore,  experimental  parameters  were  adjusted  to  produce  background  optical  densities  less 
than  this  value — ^with  an  attendant  decrease  in  dynamic  range. 

To  minimize  penumbra,  the  film-to-jet  distance  was  made  as  small  as  possible.  The  minimum  object- 
to-film  distance  (r^j  =  4.4  cm)  was  imposed  by  the  dimensions  of  the  test  chamber.  The  x-ray  tube-to-film 
distance  (r^j  =  42.4  mm)  was  a  compromise  between  (1)  penumbra  considerations,  which  are  reduced  with 
increasing  r^j  (equation  6);  and  (2)  achieving  the  desired  optical  densities,  which  fall  off  as  It  was 
found  that  jet  details  with  good  signal-to-noise  ratio  were  obtained  when  image  optical  densities  fell  in 
the  range  OD  =  0.9  to  2.  Substitution  of  these  r^j  and  r^^  values  into  equation  6  yields  a  penumbra  value 
of  0.31  mm.  Since  the  actual  resolution  is  typically  five  times  better  than  the  penumbra,  the  goal  of 
resolving  core  structure  details  to  less  than  0.1  mm  was  achieved. 
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3.  EXPERIMENTAL  SETUP  AND  OPERATION 


The  experimental  setup  and  detaUs  of  the  major  components  are  shown  in  Figures  1  and  2.  The  test 
chamber  is  a  20.3-cm-long  steel  cylinder  with  a  5.72-cm  LD.  The  chamber  contains  two  pairs  of  opposing 
ports  at  90°  to  each  other.  One  pair  of  ports  contains  windows  (9.8  cm  x  3.5  cm)  which  provide  access 
for  radiography  and  photography.  The  front  window  of  this  pair,  through  which  it  was  desired  to  have 
visible  access,  was  made  of  1.9-2.5-cm-thick  polycarbonate  (Lexan™).  Lexan™  is  visually  transparent 
and  has  a  relatively  high  strength  vs.  x-ray  absorption  cross-section  ratio.  Since  visible  access  was  not 
required  through  the  rear  port,  it  was  fitted  with  a  polypropylene  window.  (Polypropylene  is  more 
resilient  than  Lexan™.)  The  radiographic  film  and  intensifying  screens  were  encased  in  a  thin,  light-tight 
plastic  envelope  and  mounted  between  the  polypropylene  window  and  a  Lexan^"*^  flange  as  shown  in 
Figure  2.  This  assembly  was  clamped  together  by  the  port’s  mounting  bolts.  A  rubber  pad  was  placed 
between  the  Lexan™  flange  and  the  film  pack  to  uniformly  distribute  the  pressure  on  the  intensifying 
screens.  The  other  ports  included  (1)  a  clear  Lexan^“  window  (7  cm  x  1.9  cm),  through  which  the  jet 
could  be  illuminated  with  visible  light  for  standard  photography;  and  (2)  a  spark  plug,  gas  port,  and  the 
thermocouple  port  (not  shown). 

Two  methods  were  employed  for  providing  high-temperature  and  pressure  ambient  gas.  The  first 
method  (test  series  A)  involved  delivering  nitrogen  to  the  test  chamber  from  a  high-pressure  gas  reservoir 
via  a  particle  bed  heater  (Bitk  and  McQuaid  1993).  The  second  method  (test  series  B— E)  involved  igmting 
a  pressurized  stoichiometric  mixture  of  H2/Air/Ar  in  the  test  chamber  with  an  automotive-t)q)e  spark  plug. 
(Argon  was  used  in  the  mixture  to  facilitate  ignition  and  prevent  detonation.)  The  temperature  (Tc)  and 
pressure  (Pq)  of  the  resulting  combustion  products  reached  near  adiabatic  limits  in  approximately  5  ms, 
then  dropped  due  to  heat  losses  to  the  chamber  walls  (Figure  3).  By  controlling  the  injection  timing  with 
respect  to  ignition,  different  (Pq,  Tc)  values  for  the  same  ambient  gas  density  (pg)  could  be  produced. 

Chamber  (P^.)  pressure  was  measured  with  a  piezoelectric  pressure  transducer  (Kistler,  Model  601B1). 
Measurement  of  the  test  chamber  temperature  depended  on  the  method  used  to  generate  the  ambient  gas 
conditions.  In  test  series  A,  where  the  chamber  was  loaded  with  nitrogen  through  the  particle  bed  heater, 
the  gas  temperature  (T^)  was  recorded  using  an  exposed  0.05-mm  chromel-alumel  thermocouple.  In  test 
series  B-E,  where  the  stoichiometric  H2/Air/Ar  mixture  was  ignited,  T^.  was  calculated  from  the  measured 
Pc  at  the  time  of  the  x-ray  pulse.  The  calculation  is  based  on  computing  the  adiabatic  pressure  (Pca) 
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2Ho+1Alr+3Ar 


Legend: 

PT  -Ki8tler601B1  Pressure  Transducers 
LW  -  Clear  Lexan  VVirKJow  (25nirT)  Thick) 
P  -  12.7mm  Polypropylene 
IS  -  3M  Trimax  8  Intensifying  Screens 
F  -Dupont  NDT-57  Rim 
R  -Rubber  Pad 
L  -20inm  Laxan  Block 


Figure  2.  Test  chamber  detail. 
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and  temperature  (Tq^)  for  the  mixture  using  the  Blake  thennodynamics  code  (Freedman  1982).  Since  the 
chamber  volume  is  constant,  and  approximately  9  mol  of  gas  products  are  generated  for  every  10  mol  of 
reactant,  T^-,  can  be  calculated  from  the  initial  gas  temperature  (Tj^q)  ^nd  pressure  (Pco)>  ^  modified 
equation  of  state  for  the  gas,  viz., 

Tc  =  (10/9)  X  [(Pc  X  Tco)/Pco]  x  [l  +  K  (Pc-Pco)/(Pca-Pco)1  >  (7) 

where  K  is  a  correction  factor  that  scales  Tc  over  the  range  Tco  to  Tc^  for  Pc  between  Pco  and  Pca- 
(K  =  0.058,  0.057,  0.055  for  [Pcq.  PcaI  =  tl-379.  10.02];  [1.724,  12.56];  [2.069,  15.12]  MPa  and  [Tcq, 
Tca]  =  [293, 2503];  [293, 2507];  [293,  2509]  K,  respectively.)  Typically,  the  maximum  Pc  obtained  was 
90%ofPcA- 

Methyl  iodide  was  injected  vertically  upward  into  the  test  chamber  via  an  injector  in  the  base  of  the 
test  chamber.  The  injector  is  pneumatically  actuated  and  utilizes  hydraulic-pressure  amplification  to  drive 
the  liquid  (Biik  and  Reeves  1987).  The  pressure  exerted  on  the  liquid  during  the  injection  was  measured 
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with  a  piezoelectric  pressure  transducer  (Kistler,  Model  601  Bl).  A  photo-interrupter  (based  on  a 
GE  H11A2  photodiode)  was  arranged  so  that  a  comb-like  mask  attached  to  the  injector’s  pushrod  passed 
through  it.  The  (square  wave)  signal  from  this  configuration  was  used  to  monitor  the  injection  velocity 
(Figure  3)  and  trigger  the  flash  x-ray  unit.  An  output  signal  associated  with  the  triggering  of  the  x-ray 
tube  was  recorded  in  order  to  ascertain  the  injection  parameters  at  the  time  the  x-ray  image  was  taken. 

In  addition  to  the  flash  x-ray  photograjAs,  corresponding  images  based  on  visible  light  scattering  were 
obtained  with  an  intensified-CCD  camera  in  test  series  A.  (In  test  series  C— E,  the  intense  radiation  from 
the  hot  gas  products  washed  out  the  camera  images.)  The  camera  employed  (EEV)  has  a  578  x  244-pixel 
array  and  was  equipf)ed  with  a  TV  zoom  lens  (Fujinon,  Model  C6xl7.5B).  To  obtain  images  that  froze 
the  fluid  motion,  the  f-stop  of  the  lens  was  set  to  1.8,  and  an  exposure  setting  of  1  ps  was  chosen.  The 
camera’s  intensifier  gain  was  set  to  maximize  the  image’s  intensity  range  at  these  settings.  The  gain 
setting  (approximately  2.5  out  of  10)  was  established  by  trial  and  error  during  preliminary  testing.  (Note 
that  exposure  is  an  exponentially  increasing  ftinction  of  the  gain.) 


To  record  the  same  field  of  view  as  the  x-ray  image,  the  camera  viewed  the  injection  via  a  3-mm-thick 
mirror  (aluminum-coated  glass)  in  the  path  of  the  x-ray  beam  and  at  45°  to  its  axis.  Preliminary  testing 
showed  that  the  minor  did  not  degrade  the  x-ray  image.  The  camera  was  also  rotated  90°  to  its  normal 
(horizontal)  orientation  to  take  best  advantage  of  the  CCD  element  aspect  ratio.  The  effective  resolution 
of  this  configuration  was  approximately  0.2  mm^ixel.  The  camera  was  triggered  simultaneously  with  the 
flash  x-ray  pulse,  and  an  output  corresponding  to  exposure  onset  and  duration  was  recorded.  Thus,  the 
jet  core  (recorded  in  the  x-ray  image)  and  the  spray  envelope  (recorded  via  the  camera)  can  be  superposed 
to  establish  the  relationship  between  the  two. 

4.  IMAGE  PROCESSING 

The  x-ray  film  exposed  during  an  experiment  was  developed  with  an  automatic  processor  (DuPont, 
Model  NDT  100)  at  normal  settings.  The  jet  appears  as  light  shades  of  bluish-gray  on  a  dark-gray 
background  in  the  resulting  transparencies.  The  dynamic  range  of  the  optical  densities  in  the 
transparencies  typically  ranged  from  0.5  to  2.5.  This  corresponds  to  a  (digital)  gray-level  range  of 
about  100.  Digital  images  of  these  transparencies  were  obtained  with  a  PC-based  optical  scanner 
(Artiscan,  Model  8000C).  The  transparencies  were  well  reproduced  at  a  scanner  resolution  of 
200  pixels/cm.  Scan  control  and  image  processing  was  done  with  Aldus  Photostyler.  All  digitized 
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images  were  handled  in  an  8-bit  gray  level  format  (TIF).  To  maximize  the  dynamic  range  in  the  digital 
recording,  the  scanner’s  exposure  was  set  to  yield  a  maximum  gray  level  value  for  the  core  adjacent  to 
the  injector  while  maintaining  background  values  near  0  ("black").  These  images  were  then  processed  to 
yield  images  with  256  gray  levels  (0-255).  This  pseudoexpansion  of  the  dynamic  range  facilitated  image 
analysis  and  reproduction. 

In  the  CCD  camera  images,  the  jet  also  appears  as  light  shades  of  gray  on  a  dark-gray  background. 
The  digital  image  created  by  the  camera  was  converted  to  a  video  signal  and  recorded  on  videotape.  It 
was  then  digitized  with  a  PC-based  frame-grabber.  Direct  transfer  of  the  digital  image  to  a  computer  is 
possible,  but  this  procedure  proved  unreliable.  The  difficulty  was  presumably  due  to  the  noisy 
electromagnetic  environment  created  by  the  experiment 

5.  TEST  MATRIX 

The  test  parameters  are  summarized  in  Table  2.  The  present  work  is  restricted  to  the  study  of  effects 
of  variation  in  ambient  gas  temperatures,  pressures,  and  densities.  The  armular  jets  had  smaller  thicknesses 
and  velocities  than  the  full  cone  jets  due  to  experimental  constraints.  In  all  tests,  the  radiograph  was  taken 
after  the  injection  pressure  reached  steady  state.  Because  only  one  radiograph  was  obtained  per  test,  a 
certain  redundancy  was  incorporated  into  the  test  matrix  to  statistically  validate  the  results.  A  limited 
number  of  tests  with  0.5-mm  and  2-mm  jets  were  also  conducted.  These  tests  provided  a  measure  of 
liquid  density  distribution  vs.  gray  level.  They  also  indicated  that  the  experimental  trends  observed  with 
the  1-mm  jets  span  this  jet  diameter  range.  As  expected,  the  core’s  penetration  distance  scaled  with  jet 
thickness. 

6.  RESULTS 

6. 1  General  Observations.  Figures  4, 5,  and  6  represent  typical  flash  x-ray  and  CCD  camera  images 
obtained  in  the  experiments.  They  demonstrate  the  success  of  radiography  in  revealing  the  liquid  core 
structure  under  all  ambient  conditions.  Furthermore,  they  show  that  the  core  is  not  a  well-defined 
stracture  whose  diameter  linearly  diminishes  with  downstream  distance.  Rather  it  is  a  zone  in  the  jet 
whose  liquid  or  liquid/vapor  concentration  per  unit  volume  is  much  higher  than  that  of  the  overall  spray. 
Many  of  the  experiments  show  accumulation  of  CH3I  liquid/vapor  far  downstream.  This  is  expected  since 
it  is  a  closed  chamber.  We  estimate  that  an  x-ray  absorption  path  length  through  0.1 -mm  CH3I  in  liquid 
form  is  equivalent  to  a  path  length  of  3  to  10  mm  in  the  accumulated  gas. 
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Table  2.  Test  Parameters®’  *’ 


Series 
(no.  of 
tests®) 

Pg, 

(g/cm^) 

Pc 

(MPa) 

Tc 

(K) 

Pl/Pg 

Po 

(m/s) 

Rec  X  10"'* 

Lc 

(cm)  (calc.) 

A’’ 

IH 

(3,2) 

6.48 

293 

55-94 

19.2-32.8 

6.23 

(8,2) 

13.79 

623 

55-89 

12.8-20.7 

6.23 

■■ 

Mi 

(2,2) 

0.0141 

mBm 

293 

ISS 

74 

4.9 

14.0 

0.0141 

293 

wm 

95 

7.3 

14.0 

(11,1) 

0.0141 

4.9-9.0 

1,187-2,225 

161.6 

75-87 

1. 3-2.3 

14.0 

D® 

Ml 

■■■ 

(3,0) 

8.0 

1,561 

75 

2.0 

10.28 

2,028 

88 

2.2 

0.0176 

11.45 

2,273 

129.5 

77 

1.5 

12.5 

M 

(2.0) 

■aj2UH 

12.1 

1,978 

108 

85 

2.2 

11.4 

mmm 

12.88 

2,119 

108 

77 

1.9 

11.4 

^  (circular  jet,  annular  jet). 

^  Injection  into  N2  (uq  =  38-44  m/s  for  annular  jets). 

^  Injection  into  2H2  +  IO2  +  4N2  +  3Ar. 

^  Injection  into  the  combustion  {»“oducts  (2H2O  +  4N2  +  3Ar)  of  test  series  B  gas. 

®  Same  as  C,  but  with  precombustion  pressure  of  1.724  MPa. 

^  Same  as  C,  but  with  precombustion  pressure  of  2.069  MPa. 

^  Rcq  =  Pg^^Pg  Reynolds  number  in  the  gas). 

^  Net  injection  pressures  (Pl-Pq)  ranged  from  8.2  MPa  to  26.2  MPa  for  circular  jets,  and  4.8  to  7.6  MPa  for  annular  jets. 
Discharge  coefficient  values  range  from  0.58  to  0.68.  Liquid  Reynolds  number  (Rcl  =  Pl^o^Pl)  circular  jets  range 

from  2.5  X  10^  to  3.  X  lOl 
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4a.  Nonevaporating 
Uq  =  90  m/s 
Pc  =  6.48  MPa 
Tc  =  293  K 
Pl/Pg  “ 


4b.  Evaporating 
Uq  =  68  m/s 
Pc  =  13.79  MPa 
Tc  =  623  K 
Pl/Pg  ”  ^ 


Figure  4.  Visual  and  x-rav  records  of  evaporating  and  nonevanorating  sprays  in  test  series  A. 
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5a. 

Uq  =  74  m/s 
=  1.38  MPa 
Tc  =  293  K 
Pl/Pg  =161.6 


5b. 

Uq  =  86  m/s 
Pc  =  6.48  MPa 
Tc=  1,581  K 
Pl^Pg  “  161.6 


5c. 

Uq  =  86  m/s 
Pc  =  8.97  MPa 
Tc  =  2,225  K 
Pl^Pg  ”  161.6 


5d. 

Uq  =  88  m/s 
Pc  =  10.28  MPa 
Tc  =  2,028  K 
Pl/Pg  =  129.5 


X-rav  records  of  jets  in  test  series  B-D:  (a)  nonevaporating,  (b)  subcritical  evaporation, 
(c)  transcritical,  and  (d~)  supercritical. 
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6a. 

Visual  record  of  top  half  of  6b. 
(Image  is  processed.  See  text  for 
discussion.) 


6b. 

Up  =  43  m/s 
Pc  =  6.48  MPa 
Tc  =  293  k 
Pl^Pg  ~ 


6c. 

Uq  =  38  m/s 
Pc  =  13.79  MPa 
Tc  =  623  K 
Pl^Pg  “ 


Figure  6.  Visual  and  x-rav  records  of  annular  jets. 


6d. 

Uq  =  40  m/s 
Pc  =  1.38  MPa 
Tc  =  293  K 
Pl^Pg  ”  161.6 


6e. 

Uq  =  38  m/s 
Pc  =  8.83  MPa 
Tc  =  2,189  K 
Pl^Pg  “  161.6 


Figure  6.  Visual  and  x-ray  records  of  annular  jets  (continued). 

Figures  4  and  5  show  that  within  a  few  jet  diameters,  the  cores  of  the  fuU-cone  sprays  become  helical, 
with  average  pitch  decreasing  with  travel  downstream.  This  is  also  noticeable  in  the  CCD-imaged 
photograph  of  the  evaporating  jet  (Figure  4b).  In  a  prior  test  program  (Biric,  McQuaid,  and  Bliesener 
1992),  we  observed  intermittent  helical  structure  in  the  high-speed  cinemagraphic  records  of  ethanol  jets 
injected  into  ambient  gas  at  a  temperature  and  pressure  above  the  critical  point  of  ethanol. 


The  core  structures  of  evaporating  jets  differ  greatly  from  those  of  nonevaporating  jets.  The  core 
penetration  distance  decreases  with  an  increase  in  ambient  gas  pressure  and  density  as  expected  (Table  2). 
Surprisingly,  however,  when  sprays  evaporated  at  temperature  and  pressure  above  the  critical  point  of 
CH3I,  but  at  the  same  global-ambient  gas  density,  core  penetration  increased  with  pressure  (in  the  range 
tested).  (It  should  be  noted  that  the  actual  critical  pressure  in  our  tests  will  be  higher  than  that  of  pure 
CH3I  because  the  ambient  gas,  which  dissolves  into  CH3I,  has  a  large  fraction  of  H2O  with  a  critical  point 
of  [647  K,  21.5  MPa].) 


Figure  6  shows  the  results  from  a  test  with  an  annular  configuration.  (Note:  In  Figure  6a,  only  one 
boundary  of  the  jet  was  illuminated  and  recorded.  The  lit  boundary  was  flipped  vertically  and  appended 
to  the  original  to  produce  the  image.)  The  core  emerges  from  the  injection  port  as  a  cylindrical  shell 
composed  of  streamwise  ligaments.  The  shell  tends  to  collapse  inward  after  a  short  travel  (1-2  shell 
diameters),  and  more  quickly  at  higher  ambient  pressures  and  temperatures.  When  magnified,  many 
ligaments  exhibit  the  helical  structure  observed  in  the  cores  of  the  full  cone  sprays. 
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6.2  Physical  Interpretation.  Helical  instabilities  of  the  type  observed  have  been  predicted  theoretically 
(Mattingly  and  Chang  1974)  and  seen  experimentally  in  both  single-phase  submerged  jets  and  in 
two-phase  liquid  jets  (Mattingly  and  Chang  1974;  Hoyt  and  Taylor  1977).  They  stem  fiom  dynamic 
instabilities  which  are  the  result  of  momentum  fluctuations  and  have  wavelengths  on  the  order  of  a  few 
jet  diameters.  Theories  which  have  been  used  to  explain  these  phenomena  are  based  on  linearized 
equations  of  motion  with  an  imposed  perturbation  of  the  form  exp[in(t)  +  ia  (x-ct)],  where  x  and  (j)  are  the 
longitudinal  and  the  azimuthal  cylindrical  coordinates.  When  n  =  1,  the  jet  trajectory  will  become  helical. 
If  a  is  a  complex  number,  the  perturbation  may  grow  exponentially  with  distance  (i.e.,  spatial  instabilities). 
The  perturbations  oscillate  in  both  time  and  space.  Our  experimental  results  can  be  explained  in  terms 
of  such  spatial  instability.  Furthermore,  the  theories  indicate  that  the  helical  instabilities  (n  =  1)  dominate 
symmetric  instabilities  (n  =  0)  when  the  jet  profile  has  more  Gaussian  than  "top  hat"  character.  This  fits 
the  observations  that  at  higher  ambient  gas  densities  (test  series  A)  the  core  becomes  helical  at  shorter 
distances  than  at  lower  densities  (test  series  B).  For  the  same  global  gas  density,  the  core  transitions  to 
a  helical  structure  more  rapidly  when  the  ambient  gas  is  hotter.  This  is  probably  due  to  higher  local  gas 
densities  near  the  jet  origin  under  these  conditions.  The  helical  trajectory  pitch  decreases  with  travel  as 
the  velocity  decreases  due  to  liquid  momentum  transfer  to  the  ambient  gas. 

While  helical  instabilities  develop  in  the  core  a  few  diameters  from  the  origin,  capillary  instabilities 
associated  with  surface  tension  are  responsible  for  the  conically  shaped  spray  envelope.  These 
perturbations  are  axisymmetric  (Reitz  and  Bracco  1986;  Lin  and  Kang  1987)  and  dominate  near  the  core 
origin.  Such  instabilities  have  wavelengths  which  are  much  smaller  than  the  jet  diameter.  If,  as  a  result, 
liquid  droplets  strip  off  the  jet  periphery,  the  liquid  core  is  likely  to  decrease  in  diameter  with  distance 
from  the  origin  until  it  becomes  helical.  Indeed,  such  behavior  was  observed  in  all  tests  except  those 
conducted  with  ambient  conditions  above  the  critical  point  (series  D  and  E). 

The  increase  in  the  penetration  distance  with  increase  in  pressure  and  temperature  above  the  critical 
point  can  be  explained  if  it  is  assumed  that  the  core  surface  (including  the  surfaces  of  liquid  blobs  within 
the  core)  achieve  critical  temperatures  within  a  short  distance  from  the  origin.  In  this  case,  the  surface 
tension  vanishes  and  the  Ohnsorge  number  becomes  very  large.  Indeed,  theoretical  predictions  indicate 
that  when  the  Ohnsorge  number  is  large  atomization  is  slow.  Atomization  is  a  relatively  effective  pathway 
for  the  transfer  of  momentum  from  the  core  to  the  ambient  gas,  and  it  increases  the  liquid  surface  area, 
facilitating  evaporation.  Slower  atomization  could  produce  the  unexpectedly  long  axial  core  penetration 
observed  in  test  series  D  and  E. 
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The  inward  collapse  of  the  annular  jet  core  is  presumably  due  to  a  low-pressure  zone  being  established 
inside  the  annulus.  The  entrainment  of  the  ambient  gas  by  the  inner-atomizing  surface  of  the  core  could 
produce  such  a  low-pressure  zone.  At  higher  gas  densities,  the  entrainment  is  more  effective  and  results 
in  faster  collapse.  In  high-temperature  tests,  it  is  probable  that  the  gas  inside  the  core  has  a  lower 
temperature  (and  higher  density)  than  the  global  environment.  This  would  explain  the  jet  collapse  despite 
the  low-global  gas  density  (see,  for  example.  Figure  6e.) 

While  the  foregoing  discussion  provides  interpretations  for  the  observed  trends,  the  scope  of  the 
analysis  is  limited.  Other  physical  factors,  such  as  the  effect  of  the  gas  viscosity  (which  becomes  higher 
with  temperature)  or  additional  processes  related  to  thermal,  mass,  and  momentum  transfer  at  high 
pressures,  may  also  be  important. 

7.  PRACTICAL  IMPLICATIONS 

The  experimental  results  show  that  the  liquid  core  structure  is  affected  by  heat  transfer  at  hi^ 
pressures.  The  penetration  of  the  core  depends  on  ambient  conditions  in  a  manner  not  accounted  for  by 
present  theories.  Therefore,  spray-combustion  codes  may  give  erroneous  results.  The  experimental  results 
indicate  that  the  theories  may  underestimate  the  spray  penetration  at  high  pressures  and  temperatures. 
These  observations  have  particular  relevance  to  the  high-intensity  spray-combustion  process  in  guns. 
Because  core  structures  are  inherently  unstable  and  fluctuate  with  time,  the  fluctuations  can  trigger 
combustion  noise  that  couples  with  acoustic  modes  in  the  chamber  and  result  in  pressure  fluctuations.  The 
situation  may  be  aggravated  if  the  pressure  excursions  periodically  transit  through  critical  levels.  Pressure 
imbalances  on  the  core  of  annular  jets  in  guns  may  also  result  in  combustion  instabilities.  However,  this 
possibility  requires  further  investigation.  Gun-injector  geometries  are  much  more  complex,  with  much 
thicker,  higher-velocity  jets  than  those  observed  in  the  present  experimental  study. 

8.  CONCLUSIONS 

•  Flash  x-ray  radiography  has  been  successfully  employed  to  resolve  small  (0. 1  mm)  details  in  the 
cores  of  evaporating  CH3I  jets  at  elevated  pressures  and  temperatures. 

•  It  was  observed  that  the  core  structure  is  strongly  affected  by  heat  transfer. 
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•  The  core  stracture  is  inherently  unstable,  fonning  helical  structures  whose  wavelengths  are  larger 
than  the  jet  diameter. 

•  The  penetration  of  all  cores  decrease  with  ambient  gas  density.  For  the  same  gas  density,  the 
lengths  of  evaporating  cores  decrease  with  pressure,  reach  a  minimum  at  some  critical  pressure,  then 
increase  with  further  rise  in  pressure.  The  minimum  corresponds  to  a  transition  to  a  temperature 
and  pressure  regime  above  the  themodynamic  critical  point  of  the  liquid. 

•  The  results  indicate  that  core  fluctuations  are  a  potential  source  of  instability  in  combusting  sprays. 


Intentionally  left  blank. 
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